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МАТНЕМАТІСАЇ, МОРЕШІХС ОЕ ТЕМРЕКАТОВЕ ЕТКІ.О5 ІХ СКО55- 
ВОКРЕК АВЕКАЗ АТ ЕГЕСТКОКІС КАРІАТІОК УУКІ.РІХС 


У статті вперше побудована математична модель розрахунку полів температури в довільних областях 
при електронно-променевому зварюванні у вигляді крайової задачі математичної фізики для пароболічного 
рівняння теплопровідності з граничними умовами Неймана. Тепловий потік у тілі при зварюванні моделювався 
точковим джерелом тепла, що рухається по контуру тіла з сталою швидкостью за допомогою функції Дірака. 
Було побудоване нове інтегральне перетворення для двовимірного кінцевого простору, із застосуванням якого, 
а також методів кінцевих елементів і Гальоркіна, знайдено температурне поле у вигляді збіжного ряду. 
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Вступ 

У виробництві зварних конструкцій 
все більшого поширення набувають туго- 
плавкі, жароміцні, антикорозійні 1 радіа- 
ційно-стійкі матеріали, для зварювання 
яких потрібні особливі методики: такі як 
електронно-променеве зварювання, при 
якому температура активної робочої зони 
досягає в тисячу разів більших показників, 
ніж при традиційних способах. 

Даний вид зварювання широко засто- 
совується при зварюванні високоміцних 
легованих сталей і сплавів на титановій 
основі, а також таких металів як молібден, 
тантал, ніобій, вольфрам, цирконій, бери- 
лій при точній обробці і зварюванні різних 
мікродеталей. Використовується у таких 
галузях як ракетобудування, ядерна енер- 
гетика, точне приладобудування, мікро- 
електроніка та багатьох інших. 
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Постановка проблеми 

Залежно від конструкційних особли- 
востей установки, електронно-променеве 
зварювання може проводитися  перемі- 
щенням променю щодо зафіксованої дета- 
лі. Тепловий потік у виробі при зварюванні 
моделюють зосередженими,  розподіле- 
ними по поверхні або об'ємними джере- 
лами тепла. 

Труднощі вивчення фізичних явищ у 
зоні впливу електронного променю обхо- 
дять шляхом введення певного джерела 
теплоти 1 використання теорії теплопровід- 
ності. Такі підходи в ряді випадків дають 
можливість швидше отримати методики 
інженерних розрахунків процесу, ніж 
детальний аналіз фізичних явищ. Тому 
подальший розвиток методів розрахунку 
температурних полів зварних виробів є 
актуальною  науково-технічною  пробле- 
мою, що дозволить отримати більш ефек- 
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тивні 1 обгрунтовані рішення про напру- 
жено-деформований стан при зварюванні. 

Ось чому до числа проблем, що 
представляють великий теоретичний |і 
практичний інтерес, належить проблема 
вивчення температурного поля у тілах при 
електронно-променевому зварюванні. 

У даній роботі будемо розглядати 
електронно-променеве | зварювання як 
спосіб з'єднання тіл довільної конфігурації. 

Аналіз останніх досліджень |і 
публікацій 

Зарубіжними дослідниками для моде- 
лювання теплового потоку при зварюванні і 
розрахунку температурних полів звареної 
конструкції застосовується модель об'ємного 
джерела тепла |3|, у якій тепловий потік, 
розподілений згідно з законом Гаусса, діє 
всередині подвійного еліпсоїда. Однак засто- 
сування цієї моделі вимагає від дослідника 
попереднього знання характерних розмірів 
зварювальної ванни (довжини, ширини, гли- 
бини), що пов'язано з необхідністю прове- 
дення натурних експериментів для обраних 
режимів зварювання. 

Для визначення температурного поля 
при зварюванні широко використовують 
як аналітичні (1-41, так і чисельні методи 
15-6). Однак, складність термомеханічних 
процесів, що протікають при зварюванні в 
малих за розміром зонах впливу лазерного 
променю на матеріал, часто унеможлив- 
лює розрахунок температурних полів для 
конструкції заданої геометрії. 

Подальший розвиток методів розра- 
хунку температурних полів зварних виробів є 
актуальною науково-технічною проблемою, 
що дозволить отримати більш ефективні і 
обгрунтовані рішення про напружено-дефор- 
мований стан при зварю-ванні. 

Тому в даній роботі пропонується 
метод розрахунку температурних полів у 
довільних областях при  електронно- 
променевому зварюванні. 

Мета статті 

Побудова нової математичної моделі 
розрахунку температурних полів у 
довільних областях при електронно-проме- 
невому зварюванні у вигляді крайової 
задачі математичної фізики, а також 
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знаходження рішень отриманої крайової 
задачі. 

Викладення основного матеріалу 
дослідження 

Розглянемо розрахунок температур- 
ного поля при з'єднанні тіл довільної 
конфігурації електронно-променевим зва- 
рюванням. Нехай Дє К'- обмежена 


область із замкненим  кусково-гладким 
контуром Г, а 7 зовнішня одинична 
нормаль до Г, то тепловий потік у тілі при 
зварюванні будемо моделювати точковим 
джерелом тепла, що рухається по контуру 
Г зі швидкістю У з інтенсивністю 0. 

Теплофізичні властивості тіла не 
залежать від температури. У початковий 
момент часу температура циліндра 
постійна Су, а на зовнішній поверхні відо- 
мі значення теплового потоку С, (ж у) 

З урахуванням прийнятих допущень, 
математично задача визначення темпера- 
турного поля Т(х, ум) складається в 
інтегруванні диференціального рівняння 
теплопровідності в області П: 


0-д |аг-у: 
00 078 070 ОР 

за| 299-953 
ОЇ дх'" ду сріТ.х - бо) 


тах 


з початковою умовою 


д(х, у, 0) ш 0, (2) 
і граничною умовою 

00 

72) «лю 

дя |, у ) З) 


де цЕЗ у,1) - температура 
тіла; у - щільність матеріалу; с - питома 
теплоємність; 


А 
2. - коефіцієнт теплопровідності; а - - и 
коефіщієнт температуропровідності; 
9(х а: дельта функції 
Та У, г) б Со 


тоб. 


тах 


Дірака; 


- відносна температура 
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тіла; Ах, у)є ка - криволіній- 
ний інтеграл; О - початкова точка знаход- 
ження точкового джерела тепла; Р - по- 
точна точка знаходження точкового 
джерела тепла. 

Для розв'язання крайової задачі (1)- 
(3) застосовуємо інтегральне перетворення: 


Га) | ебо ма )УМо 0 


де фі, у, дк. ш. - власні функції і власні 
значення. 

Класична проблема власних значень і 
власних функцій формулюється як задача 
про визначення значень числових парамет- 
рів (власні значення) //, і функцій (власні 


функції) фі, З о які тотожно нерівні 
нулю в області |) та задовольняють 


рівнянню: 
д'ф д/ф 
фозУЗЦоФУЙ 
а бу оф (5) 
і граничний умові: 
СИ 
п |, (6) 


де фо Узі, ) сС(Фуз Їм, у) є С(Р): 


дах, У) є С(Р), Ма, | а |є2 У; 


дих 
оон кь пара, за, -муль 
а ; | а аз , 
Ох "ОЖо 
тиїіндекс, компоненти якого Є цілі 
невід'ємні числа. 
Знайдемо власні значення щ і власні 


функції фі, У, и,) із розв'язку задачі (5)- 
(6) за допомогою методів кінцевих 
елементів і Гальоркіна. Для цього зробимо 
розбиття області на симплекс елементи 


(рис. 1): 
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1 (ха, мі) М 
Х п 


Рис. 1. Трикутний елемент порядку. 


ф.(х, у) всередині 
елемента виражається через 


Тоді функція 
симплекс 
функції форми М,, МУ, и М, із відомими 
значеннями  ф, ф,, Ф; у вершинах 


трикутника: 


ф.(х, у) я Мф -- М,ф, - Мф, - 
м в (т 


де ІМ ДАТМ, М МОЇ з фр фу 
нижній індекс (е) означає довільний 
симплекс елемент. 

Для іго вузла (і - 1, 2, 3) функції 
форми мають вигляд: 


М (ху) є за чу нсу), 


де 
1 
5, зв з уд ох) З (ху зх ОУ» з У 


а ЗКУ ХУ Б, ЗУ) Укр ХХ 


5, - площа трикутного елемента; (х|, У|), 


(х»,У2) 1 (хз Уз) - координати його 
вершин; і, /,К - послідовна нумерація 
вузлів симплекс елемента при обході їх 
проти годинникової стрілки. 

Для визначення функцій форми М,, 


М,, МУ, до зручно використовувати І, 


і 
координати всередині симплекс елемента, 
які визначаються відношенням площі 
трикутника, | утвореного | точкою 1 
стороною, протилежною вершині і, до 


загальної площі трикутника: 
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1 

І "зу на -вуюм о 39 5 
1 

І, "зу НОРМ з 
1 

І, 255 Таб чу У р ; 


Де Ср Уа Как 1973/1783 091 392 
у ЗХ37Х) У ЗХ7Х У Хо ЗД: 

У разі лінійного симплекс елемента, що 
містить отри вершини, функції форми 
збігаються з відповідними І, координатами: 

М, «5, МУ» 2 5, Мз є з. 

Для представлення задачі на власні 
значення в матричній формі використаємо 
метод Гальоркіна. Підставимо в рівняння 
(5) наближений розв'язок (7), тоді 
отримаємо рівняння: 


ЗЕ т 
ооо || М Че, 
ре А МУК; 
и. ЇМ. (фа 0 (9) 
Множення у лівій частині рівняння 


(3) на функцію форми |М,|, та інтегру- 
вання по елементу е дає: 


ІТ, 440) (9) 
«|| на 23) ІМ. хау 9); 


І, зм, ДЦ ОМЛІМ, 1 ахау фі). 
о. 


Інтегруючи І, пох 1у, отримаємо: 
5ов-Ї,, 
де 


ом 
ера м дучн 


Па ун їх 
бу 
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сова ура Пе 
ду ду 


Враховуючи тотожне співвідношення 


дФ 9Ф а 
є РО 


отримуємо 


мг 
- , кін 


де д/дп - похідна по зовнішньої нормалі; 
АГ - криволінійний інтеграл по межі. 
Тоді (9) набуває вигляду: 


ка І аг) - 


Оп 
Р 101М,17 «багна 
ду || бу 
ахау -(ф.)ч 


щу МЛМ ахауі9,) о (0) (10) 


Сумування по всіх елементах у (10) 


дає: 
5 Пд |З МЕ кім із 
Оп 


т. УДИМЛГМ," ахау (9) - 
Зі 


р 
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ея (ам -(0). (11) 


Помноживши перший доданок (11) 
б 
на вираз (9, Ф,) , одержуємо: 
го 


із Уме рові зона 
ЗАМ), є 
Оп 


5 Годїм.д): 


а « 


При цьому криволінійні інтеграли по 


межах Г,елементів усередині області ДР 


скорочуються, оскільки напрями обходу 
межі при обчисленні криволінійних 
інтегралів протилежні для кожної пари 
сусідніх елементів. У результаті залиша- 
ється тільки криволінійний інтеграл по 
межі Г всієї області ЇЙ. 


Тоді, враховуючи граничну умову 
Неймана, можна знехтувати | першим 
доданком у (11). Отже, (11) набуває вигляду: 


(КІ м МДФ) е (0) (2 


ку 191м, 


оїм1дТМ, мах 
ду ду 


ім3- |, М.М, ахау; 
1-31 


в: 


де 


Таким | чином власні | функції 
фо, у, УРА і власні значення //, знахо- 


дяться із (12), а формула оберненого 
перетворення має вигляд |7Ї: 
ан | р, 2, и,) 


"ру гри Г(щ.)- (3) 
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Застосовуємо до диференціального 
рівняння (1) інтегральне перетворення (4). 
У результаті отримуємо звичайне диферен- 
ціальне рівняння 


48 п 
Р - О(/) 19 


з початковою умовою 


9(0)-0. 15) 
де 


Ом) в 
: мо ср(Тих -С,) 


Шебсо зум) да пре 


ідо у м, А(х УМ 


Г-додатно орієнтований контур. 


Розв'язавши диференціальне рівняння 
(4) з початковою умовою (15), одержуємо: 


д(шуг)є | О(о м, ) ей о: 


о 


Таким чином, з урахуванням формули 
оберненого перетворення (13) отримуємо 
шукане температурне поле д(х, у, г): 


в(х, умі) - -Ув (шу»г) 


ка 


Со, у, м.) 
иа уки 


я 


Висновки 

У статті вперше побудована матема- 
тична модель розрахунку полів темпера- 
тури в довільних областях при електронно- 
променевому зварюванні у вигляді крайо- 
вої задачі математичної фізики для пара- 
болічного рівняння теплопровідності з 
граничними умовами Неймана. Тепловий 
потік у тілі при зварюванні моделювався 
точковим джерелом тепла, що рухається 
по контуру тіла з сталою швидкістю за 
допомогою функції Дірака. Було побудо- 
ване нове інтегральне перетворення для 
двовимірного кінцевого простору, із засто- 
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суванням якого, а також методів кінцевих 
елементів і Гальоркіна, знайдено темпе- 
ратурне поле у вигляді збіжного ряду. 
Знайдено розв'язок крайової задачі тепло- 
обміну | в довільних областях при 
електронно-променевому зварюванні може 
знайти застосування при моделюванні 
температурних полів, які виникають у 
тілах довільної конфігурації. 
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КЕЗОМЕ 

М.С. Вегдапук 

Маетагбісаї тодсеіїпо ої 
(етрегаїиге еїйя іп сго55-Богдег агеаз аї 
еіесігопіс гадіайоп угеїдіпо 

ІП Фе агіїсіе Їог Ше Бг5і йте а 
таФетайса! тоде! їог саїсиіацтпя, (етрега- 
сиге Пе!45 їп агіігагу боцпаєа геріопя у а 
сіозед ріесем/ізе 5пооїБ сопіойг ипаєг 
сіесітоп-Беапі  умеіЧіпє їп Ше Сагіебіап 
соогаїпаге 5у5кет їп (Бе Їогпа ої а боцпаагу 
уаїше ргобіепті Їог а рагабоїйс едиацоп ої 
Фегтаї сопдисйутну у іпіца! апа Бойп- 
Чдагу сопдаїйоп5 15 соп5(гисіед, ргоуїдеад (Паї 
Фе Фегторбузіса! ргорегіїе5 ої (Бе Боду аге 
сопягапі. Ді ап іпіца! йте, Ше іетрегаїшге ої 
Ше Боду 15 сопзкапі, апа оп Ше оціег 5игіасе 
ої Фе Боду аге Кпом/тп уаїез ої (Бе Беаї Пих 
мурісП аге сопіпиоцз5 соогаїпаге апсПоп5. 

Го (бі5 рарег, Фе реас их іп Фе Боду 
дигіпє ууеіЧіпє ма5 подеісд Бу Ше роїпі 
5ойгсе ої реа: тоуіпеє аїопе Бе Бодуз8 
сопіоиг ул а сопяапі уеіосігу апа а Кпом/п 
іпеп5їсу ц5іпеє Фе Дігас Гипсіоп. 

То 580Іуе Ше Боппдагу уаїе ргобіет, а 
пему іпіеста! (гап5їогт їог а пмо-Фітепзіопаї 
пе 5расе уга5 соп5ігисіва. Омуп уашез апа 
Феїг омуп Гипсйопя Їог (Бе іп(еєта! (гапяїог- 
тацоп Кегпе! аге Їоцпа ц5іпє Ппіе еіетепі 
апд Саїіегкіп пеоайз. Пі і8 са5е, Ше аїмі- 
зіоп ої Ше агеа іпіо 5ітріех еіетепіє м/ає 
таде. Три58, Фе ргобіеті ої Ппаїпє еідеп- 
уаїшез апа еісепіипсіоп5 у а8 гедисей їо (Бе 
аїєебгаїс ргобіета ої Ппдіпє еібепуаїшез апа 
еібепіипспопз. 

Айег арріїсацоп (о Фе обіаїед Бойпаагу 
уаїше ргобіет ої (ре сопя(гисівд пему іпіерта! 
сап5їогтацоп, ме обгаїпед Фе Сапспу ргобіет 
ууУпо5е 501шійоп5 ууеге Їоцпа апаїуцсаПу. ТПе 
обгаїпед 501ийоп ої (ре боппдагу уаїце ргобіет 
і5 (млсе сопіпиоц5іу дїНегепйакед Бу врапа! 
соогдїпаке5 апа опсе рег йте. 

Аз а гезиї, Ше (еппрегаїшге Неї іп 
агбійгагу боцпаєа геріопя млі а сіозей ріесем/і5е 
зтооі сопіоиг аї еіесітоп-Беаг ууеідтпе, 15 Їоцпа 
їп Фе Гог ОЇ сопуегеепі 5егіе5. 

Тре (етрегашге Пе! їоцпа їп агіігагу 
рБоцйпдед гебіоп5 ур а сіо5ед ріесемуїзе 
зтооіб сопіоиг аї еіесігоп-Беат ууеіФіпеє сап 
Бе иц5ед їог подиіайте (Пе (еппрегаиге Пе!45 
аг агі5е ууПеп муеіЧіпє а соп5ігисйоп ої а 
5їмеп геотеїгу. 
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